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ABSTRACT 


X-ray and ultra-violet line emission from hot, optically thin material 
forming coronal active regions on the sun may be described in terms of an 
emission measure distribution function, 0(T) ~ n e ^ dV/dT. We develop the 
relationship between line flux and 0(T), a theory which assumes that the 
electron density is a single-valued function of temperature. The sources 
of error involved in deriving 0(T) from a set of line fluxes are examined 
in some detail. These include errors in atomic data (collisionai excitation 
rates, assessment of other mechanisms for populating excited states of 
transitions, element abundances, ion concentrations , oscillator strengths) 
and errors in observed line fluxes arising from, e.g., poorly - known in- 
strumental responses. We discuss two previous analyses in which 0(T) for 
a non-flaring active region is derived, A least-squares method of Batstone 
et al . uses x-ray data of low statistical significance, a fact which 
appears to influence their results considerably. The work of Chambe, in 
which a functional form for 0(T) is deduced, uses only approximate atomic 
parameters. Replacing these with improved data, we calculate line fluxes 
using Chambe f s 0(T) and compare these with observed values. The comparison 
indicates that Chambe’s 0(T) describes the observed line emission very well. 

Two methods for finding 0(T) ab initio are developed. In the first, 
0(T) is expanded in the form 0(T) = 2 a n T n , the a n s being evaluated by 
least-squares analysis of fluxes of a number of lines. In the second, 
the integral in the expression for line flux in terms of 0(T) is replaced 
by a Gauss -La guerre quadrature sum. The coefficients are again evaluated 
by least squares. These two methods should have application not only to 
active-region plasmas, but also to hot, flare -produced plasmas. 



1. INTRODUCTION 


It is now becoming clear that the great variety of x-ray and ultra- 
violet line emission observed from coronal active regions cannot be ex- 
plained in terms of the emitting region's having a single temperature. 

This was the assumption in earlier work using data from rocket-borne instru- 
mentation (e.g. Evans and Pounds, 1968). Often, the justification for 
this was that the observed lines were emitted by a set of ions which were 
known to be produced over a fairly narrow region of temperature, so' that 
temperatures and emission measures could be deduced simply from the in- 
tensity ratio of any two lines. With the great improvement in data quality 
being achieved at the present time, this simple model is not adequate. 

It is necessary to consider models for coronal active regions in which 
there is a continuous distribution of emitting material with temperature. 

A number of recent works (Batstone et al . , 1970; Chambe , 1971; Acton 
et al . . 1972; Walker et al . , 1974) have investigated such models using 
soft x-ray observations of active regions. The x-ray line fluxes are 
analyzed to give what Batstone et al . have called "differential emission 
measure", but as we shall see this is in fact the function 0(T) = n e ^ dV/dT 
(rig = electron density, dV an element of volume,' T the electron temperature). 
In this paper, we shall consider two of these works in some detail, parti- 
cularly from the point of view of possible sources of error, to see what 
exactly has been established about the form of the distribution function, 
0 ( 1 ). We shall then develop two ab initio methods of finding 0(1) which 


have not hitherto been used. 
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2. X-RAY LINE FORMATION IN CORONAL ACTIVE REGIONS 
We shall first develop the relation between the flux of x-ray 
line emission and the function 0(1) = n e 2 dv/dl, taking such emission to 
arise from optically thin material. 

Considering the strong resonance lines that have been used in the 
works mentioned earlier, excitation of the upper level occurs predomin- 
antly by electron collisions from the ground state. For particular 
transitions, cascading from higher levels or recombination of the next 
higher ion may contribute significantly to the upper-level populations. 

When the effects of cascades and recombinations can be neglected, the 
flux (F) observed at the earth may be written as 

F = — — f n(ion) n e C 12 dv photons cm -2 s -1 (1) 

4rrr2 V 

where r is the sun-earth distance, C\2 is the rate coefficient of colli- 
sional excitation from level 1 to 2, n(ion) the ion number density, and 
the integral extends over the emitting volume, V. For a maxwellian dis- 
tribution of electron velocities, is related to the cross-section for 

collisional excitation, Q-^j by 

Clz 9l2(E) Ee ' E/kI dE c " 3 S_1 

where m is the electron mass, E the energy of the incident electron, and 
W the threshold energy for excitation. Cross-sections are often expressed 
in the Bethe approximation, for which Q 22 rs proportional to the absorption 
oscillation strength, f 12 » and to the Gaunt factor, g. This approximation 
assumes that excitations occur mostly as a result of distant electron en- 
counters, and is valid for electron energies much larger than the threshold 
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energy, Wl Following van Regemorter (1962), who uses the Bethe approxi- 
mation for Q^ 9 > we have 

o e ~W/kT f o i 

C 12 = x 1(T 3 11 P(W/kT) cm 3 s^ 1 (2) 

W /T 

where the function P is related to the Gaunt factor, g, by 

00 

P(W/kT) = J ge ”E/ kT d(E/kT). 

W/kT 

Van Regemorter has constructed an interpolation formula which approximates 
various measurements and calculations of g for s p transitions, and tab- 
ulates the resulting values of P as a function of W/kT. These indicate 
that, for positive ions, P tends to an asymptotic value of 0*20 for W ^ kT. 
From (1) and (2), we may now express F as 

F > 6.05 x 1031 hi ? . felg«Pt? n n (lon) n(H) dv 

W n(H) n (element) v e jy 

-2 -1 

photons cm s (3) 

where n(element) is the number density of all ion stages of the element 
in question, and n(H) the number density of H atoms. For weaker lines 
arising from higher levels, the right side of (3) must be multiplied by 
an appropriate branching ratio. Taking the abundance of helium relative 
to hydrogen to be n(He)/n(H) = 0.04 (e.g. Robbins et al . , 1974), we have 

n(H) = — = 0.93 n_ . 

l+2n(He)/n(H) 
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We can now write 

F = K J G(T) n e 2 dV 
V 


where 

K - 5.60 x 10 31 £12 , 

W n(H) 


( 4 ) 


C5) 


and, following a notation of Pottasch (1964), G(T) is given by 

G(T ) - ■» . e' w/kT mim . 

n(element) 


( 6 ) 


Jefferies, Orrall, and Zirker (1972) have discussed equations describing 
line emission in optically thin material which have forms similar to (4). 
Their procedure is to express n £ dV - dN e in terms of a bivariate distri- 
bution function p(n e ,T) given by 

n e dV = dN e = N e p(n e ,T) dn e dT, 

with N e the total number of electrons in the emitting volume, V. Applying 
this to (4), we have 

00 00 

F = K N e J J G ( T > n e H(ne» T ) dn e dT - ( 7 > 

o o 

The data will not in general permit us to determine the complete functional 
form of ^(n e ,T). However 5 we can find the integral of }i(n e5 T) over all 
densities providing that - n e (T) is a single -valued function of T, i.e, 
to each T, there is a single value of n e - Evidence that this is so is 
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provided by the few estimates we have of coronal active-region densities 
from white-light measurements (Newkirk 1967). Thus (7) may be expressed 
in the form 

CO 

F 53 KN e J G(T) n e (T) $(T) dT 
o 

where 

cn 

$(T) = J* n(ti e ,T) dr^. 
o 

For the purposes of our problem, we put 

0(T) - N e n e (T) ?(T) 

so that 

CO 

F = K J G(T) 0 (T) dT (8) 

o 

Comparing this with (4), we see that (8) is equivalent to putting 

0(T) = n £ 2 (T) . 

0(T) may be considered to be the function describing the manner in which 
the active region's emitting plasma is distributed with temperature. 

3. SOURCES OF ERROR 

A great deal of atomic data is required in any derivation of 0(T), 
as can be seen from the various terms occurring in Equation (3). For 
the weaker lines in the soft x-ray region, many of the parameters may be 
very poorly known, while even for the strong resonance lines used in 
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previous 0(T) analyses significant uncertainties may also exist. We shall 
discuss some of the more important sources of error. 

(i) Rate Coefficients for Direct Collisional Excitation. 

Expression (2) with van Regemorter's (1962) values for P(W/kT) 
has been extensively used in the literature for the rate of direct colli- 
sional excitation. The resulting rate coefficients agree well with calcu- 
lations using more recent cross-sections from Coulomb-Born or similar 
approximations for many of the strong resonance lines which involve s -* p 
transitions in the soft x-ray region. In particular, the asymptotic 
value of P = 0.20 for W :> kT agrees to within + 50% with values for 
2_g2 lg _ ]_ 3 2p transitions in He - 1 ike Si, Ca, and Fe ions using Coulomb- 
Born cross-sections of Blaha (1973), and for Is - 2p transitions in hydro- 
genic ions using the Coulomb-Born-Oppenheimer cross-sections of Burgess, 
Hummer, and Tully (1970). The work of Burgess et aL , which also has 
application to certain complex ions, includes the effects of exchange of 
the incident electron and the bound electron. Application of van Regmorter's 
procedure to more complex ions and transitions other than s -* p may not 
have high accuracy, however, because P(W/kT) is expected to depend on 
both ion charge and transition type (Bely and van Regemorter, 1970) . Also, 
for many x-ray transitions, we are interested in values of the collisional 
excitation cross-section close to the threshold energy, W, since, when 
W > kT (as is usually the case), only electrons in the high-energy tail 
of the maxwellian distribution are capable of exciting the transition con- 
cerned. Unfortunately, the Coulomb-Born approximation is unsuited for 
such a condition and may thus lead to incorrect results. 
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For certain ions and transitions , the occurrence of resonances in 

the collisional excitation cross-sections, arising from the formation of 

intermediate auto-ionizing states whose lifetimes determine the width of 

the resonances, may be of some importance astrophysically (e.g., Gerjuoy, 

1966). As an illustration, Gerjuoy mentions the case of a resonance 

3 

structure at 22.4 eV in the excitation cross-section of the Is 2s 
state in neutral helium. This appears to arise from the production of 

He"' in a Is 3s^ state, which auto-ionizes to form neutral He in a 

2 

2 state. Hershkowitz and Seaton (1973) have pointed to the importance 

of resonances near the threshold energy for a C III transition visible in 
the optical spectra of quasi-stellar objects. We cannot assess the 
significance of such effects for x-ray lines from the sun since there are 
not yet any data referring to resonances in the cross-sections for the 
transitions involved . 

(ii) Population of Excited States By Cascading and Recombination . 

The contribution to excited-state populations by cascading has 
been considered only approximately, but values given by Mewe (1972) indi- 
cate that cascading will be unimportant (< 10% of the direct collisional 
rate) for ls*~ - Is 2p P transitions in He-like ions and Is - 2p transitions 
in H-like ions. The same may not be true for intercombination and for- 
bidden transitions in He-like ions (although previous analyses have not 
used these lines) or for more complex ions generally. An example is the 
Ne-like Fe XVII ion, which emits a number of strong lines between 13A and 

o 

17A. Statistical equilibrium calculations for this ion by Loulergue and 
Nussbaumer (1973) suggest that the upper levels in the 2p^ - 2p-* 3s 
transition array are populated almost entirely by cascades from 2p 3p. 
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The excited states may also be populated by processes of radiative 
and dielectronic recombination. As discussed by Tucker and Koren (1971), 
the number of radiative recombinations to any given ion must always be 
less than the number of excitations from the ground state to the first 
excited level. Thus, radiative recombination as a mechanism of populating 
excited states is never more important than direct collisional excitation, 
and in fact can usually be neglected. With dielectronic recombination, 
the dielectronic capture of an electron by an ion, X , forms a doubly 
excited state in ion : 

X**" 2 (nl) + e - X +Z ' 1 (n / A / ; n'V'), 

which is stabilized by one or more radiative transitions to the ground 
state of the ion, X +Z ' :L . If the energy of the dielectronically-captured 
electron (n" &" ) can be neglected, the stabilizing transition, 

X +Z_1 (n / J e / ; n" l" ) - X +z_1 (nX ; a" &* ) , 

may be indistinguishable from the energy of the resonance -line transition, 
X +Z (n / A / ) -> X"*" 2 (nf,) , in the ion X* 2 . A discussion by Gabriel (1972) shows 
that a contribution of < 107. is made to the intensity of the Is -1 A S - Is 2p P 
transition in He-like ions by dielectronic-recombination satellite tran- 
sitions, Is 2 1 n t," -» Is 2 n" l" (n" > 3), occurring in Li-like ions. 

An approximate treatment by Tucker and Koren indicates that the contribution 
made by dielectronic recombination to line intensities in all other ions 
is less than ~ 257* of the direct collisional rate for z ^ 25. 
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(iii) Element Abundances 

Analyses of intensities of lines emitted by more than one ele- 
ment must assume values for relative abundances. Since x-ray active re- 
gions are situated in the corona, it would seem natural to assume coronal 
abundances. However, there is a possibility that the abundances of heavy 
elements will be enhanced in regions where the temperature is greater than 
the general coronal value since, for a highly ionized gas, thermal dif- 
fusion tends to drive heavy species to higher temperatures (e.g., Seaton, 
1964; Nakada, 1969). This will occur only if mixing by turbulence is 
small (a condition which seems likely to obtain in view of the presence 
of magnetic fields in coronal active regions), and if there is a temper- 
ature gradient along the magnetic-loop configurations. According to cal- 
culations by Nakada (1970), a diffusive equilibrium will be attained for, 
e.g., 0 VIII ions in a time comparable with the lifetime of an x-ray 
active region. 

(iv) Ion Concentrations 

Ionization equilibrium in a non-flaring active region has been 

assumed in all previous analyzes to deduce ion concentrations. This 

should be justified in view of the fact that, with electron densities of 
9-1 

s 10 cm - ' 5 (e.g. Newkirk, 1967), such an equilibrium will be set up very 

quickly. (Gabriel's (1972) evidence that an active region is an ionizing 
plasma must be considered doubtful since it relies upon the observed in- 
tensity of an extremely weak spectral feature.) Some existing calculations 
of ion concentrations may have to be modified slightly to take into 
account more precisely the effects of electron density on dielectronic 
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recombination rates (Summers (1973) has done this for C, 0, and Ne ions) 
and adjustments in radiative recombination rates for complex ions (Tarter, 
1971), but the effects of each should be small. However, the dielectronic 
recombination rate may show significant departures from that given by 
Burgess's (1965) formula which has been applied extensively in ionization* 
equilibrium calculations. Shore (1969) points out that the formula con- 
siderably overestimates the rate whenever the stabilizing transition in- 
volves a change in the principal quantum number. This occurs in particular 
when a H-like ion recombines to form a He-like ion. Both these ions 
feature prominently active-region x-ray spectra. 

(iv) Oscillator Strengths 

Many of the weaker x-ray lines may have poorly known oscillator 
strengths, but the errors in the calculated f-values for the strong lines 
should now be quite small as increasingly sophisticated formulations are 
applied * 

It would be clearly desirable to assess what effect each of 
the errors listed above would have on the derived functional form of 0(T) . 
However, we often have only a vague idea of the error to be placed on 
some of the quantities. To illustrate the point, inter-agreement of 
several Coulomb-Born calculations of excitation cross-sections (see (i)) 
would not necessarily imply the cross-sections (and thus the rate of 
collisional excitation) is known to high accuracy. Each calculation could 
be appreciably in error if the effect of auto-ionizing resonances is im- 
portant. If enhancement of heavy-element abundances through thermal 
diffusion cannot be neglected (see (iii)), the methods of Walker et al . 
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and Jakimiec et al . (1973), in which abundances are adjustable in the 
analysis, may have to be tried. Shore's corrected dielectronic-recom- 
bination rates (see (iv)) could be included approximately in existing 
ionization-equilibrium calculations by making empirical adjustments (e.g. 

Tucker and Koren (1972)) to the H-like and He-like ion concentrations. 

In addition to the atomic parameters, the solution of 0(T) from (8) 
will require the observed fluxes from a set of lines whose G(T) functions 
completely cover the temperature range of interest ((1-9) x 10^K, say, 
for a non-flaring active region). The observed fluxes will depend on 
the accuracy of the instrumental calibration and may, if small, be un- 
certain through count-rate statistics. A significant difficulty in the 
intensity calibration of Bragg crystal spectrometers which have been used 
in the soft x-ray region is the calculation or measurement of the inte- 
grated reflectivity coefficient, R(\). An in-flight calibration would 
obviate the need for calculating R(\), a procedure which in any case re- 
presents an ideal situation unlikely to occur in practice (see, e.g., 

Guinier, 1963). Much existing data from satellites are uncalibrated, 
however, and we must rely on theory for the calculation of R(A.) . Alexandropoulos 
et al . (1973) have investigated the change of crystal properties (R(A) 
among them) when the crystal face is momentarily heated by strong visible- 
light illumination, but it is uncertain how their results, which indicate 
significant changes in R(A) , would refer to illumination by sunlight. 

The situation seems most critical for rocket-borne detectors without in- 
flight calibration since, with many of the instruments aboard satellites, 
the crystal face is either subjected to prolonged sunlight illumination, 
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during which time R(\) assumes an equilibrium value, or (as with the 
GSFC 0S0-5 spectrometer) is shielded by a filter designed to reduce ultra- 
violet background. 

4. 0(T) FOR AN ACTIVE REGION : THE PRESENT EVIDENCE 

Solutions for 0(T) for a non-flaring active region, as we saw in ^1, 
have been claimed by various workers, but since Acton et al . (1972) and 
Walker et al . (1974) assumed that the functional form for 0(T) was 
basically that given by Chambe (1971) , the methods of analysis tried so 
far break down to only two: one given by Chambe and the other by Batstone 

et al . (1970). 

Batstone et al . have used the fluxes of five x-ray lines in a least- 

squares procedure to solve for up to four values of 0^ given by 

Tj.+AT 

H = — f 0(T) dT 

AT 

where AT = 2 x 10 6 K and T ± = 2.5 x 10 6 K, 4.5 x 10 6 K, ... . The problem, 
which can be stated in the form 

F j = aq 0i + ^2 ^2 + a 3 ^3 + a 4 ^4 CJ = 1,2,... ,5), 

is thus only just over -determined , Three active regions were observed 
by them with an uncollimated crystal spectrometer. Photon counting rates 
for individual regions were extracted by de- convolution of the data with 
the known instrumental response function- Their procedure did not strictly 
minimize the sum of squared deviations since two constraints were first 
imposed: the solutions 0^ were forced (a) to have non-negative values, and 
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(b) to decrease monotonically with T. The first of these conditions may 
appear to be reasonable, but it might have been thought that, with ob- 
servational data of adequate statistical significance , it would be un- 
necessary actually to impose this as a constraint on the analysis ♦ A 
least-squares calculation with no such conditions ought to have only pos- 
itive solutions unless there is an insufficient number of degrees of free- 
dom or the data have a poor statistical quality. 

In view of the low photon counting rates and difficulties associated 
with the de-convolution process, I have applied an unweighted least- 
squares analysis on the data exactly as presented by Batstone et al . 
without the two conditions imposed by them. In addition, I investigated 
the effects on the solutions, 0-^ , of statistical errors in the ob- 
served fluxes, F j , which I took to be Fj /Nj/N. (Nj 5:5 photon counting 
rate). Table I lists the values of 0^, together with their standard de- 
viations, for the three active regions. These are compared with values 
given by the calculation of Batstone et al . For many solutions, we have 
the physically absurd situation of negative 0^ s, but it is clear that 
this arises from the poor photon statistics since the standard deviations 
are so large. It seems that values of 0^ for T > 3.5 x 10 K are insig- 
nificantly different from zero. Data of a much higher quality are evidently 
needed for this type of procedure. 

As already discussed, Chambe (1971) assumes a functional form for 
0 (T) which, for T > 2 x 10^K, exponentially falls with T. This is alleged 
to fit the observed line fluxes very well, but in fact it is not clearly 
stated that the form of 0(T) derived actually reproduces these fluxes to 



TABLE I 


Values of 

T - (T+AT) (K) 

Active Region. It B 

P 

Active Region 2: B 

P 

Active Region 3: B 

P 


Units of 0: 10^ 


l T+AT 


0 s= ~ J 0(T) dT for Three Active Regions Seen by Batstone et al . 


(1.5-3. 5) x 10 6 (3. 5-5. 5) x 10 6 (5. 5-7. 5) x 10 6 (7. 5-9. 5) x 10 e 


15 

18.5 + 11.7 
8 

12.2 + 7.2 

16 

43.2 + 9.0 


9 

4.2 + 10.8 

4 

- 1.3 + 6.4 
8 

- 13.6 + 7.3 


7 

14.7 + 13.7 
2 

8.7 + 8.6 
6 

25.4 + 10.4 


0 

- 4.5 + 8.4 

1 

- 2.2 + 5.4 

0 

- 9.5 + 6.7 


cm”^(10 6 K) -1 . B = values from Batstone et al . (1970) 

p = values of this paper (errors are standard deviations). 
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\ 


within the observational errors . This could have been done by comparing 

CO 

observed fluxes, F obs , with the values of F ca j^ c = K j* G(T) 0 (T) dT (see 

o 

equation (8)), 0'(T) being the derived 0(T). Chambe's analysis is based 
on twenty-seven x-ray and ultra-vilet lines. The function 0(T) is de- 
rived using atomic parameters for the lines which must in many cases be 
considered very approximate. For example , values of the collisional 
excitation rates were calculated with the function P(W/kT)(see equation 
(2)) set equal to 0.3 for all x-ray (A < 100A) lines regardless of the 
transition type. As we saw, such an approximation may not be good for 
ions which are not H-like or He -like. For the Ne-like ion, Fe XVII, Chambe 
multiplies the direct collisional excitation rate, C(2p^ - + 2p ^3s) , by 
three to take account of cascading from the 2p ^3p levels. However, the 
calculations of Loulergue and Nussbaumer (1973) indicate that the cascading 
contribution will be much more significant, and a correcting factor of at 
least twenty is required. Loulergue and Nussbaumer calculate that the 
cascading contribution is also significant for the 2p ^3d - 2p^ 1 S G 

transition at 15.45A. This was totally neglected by Chambe. Finally, 

Chambe uses Bely T s (1967) calculations of the ionization equilibrium of 
iron which do not include radiative recombination and deal only approximately 
with the dielectronic recombination rate. 

In view of the possibility that the uncertainties in the atomic para- 
meters have led to inaccuracies in the derived functional form of 0(T) , I 
have investigated how well the observed fluxes of the spectral lines 
Chambe considers are reproduced by his model. A comparison was made be- 
tween the observed line fluxes, F Qbs , and the values calculated from F ca ^ c = 
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K J G(T) 0 , (T) dT, 0 / (T) being the function 0(T) given by Chambe. The 
values of KG(T) were taken from Tucker and Keren's paper where avail- 
able, except for lines due to Fe XVII for which Loulergue and Nussbaumer s 
treatment was used. Thus, element abundances as given by Tucker and 
jlo^an were assumed, as were the values given by them for the other atomic 
parameters of equations (5) and (6). There are three lines in Chambe's 
list at ultraviolet (X > 100A) wavelengths due to the ions Fe XV and 
Fe XVI. These were considered in the same way as for the x-ray lines, 
although the ultraviolet lines are marginally optically thick in their 
Doppler cores. A rough calculation, in which the extent of a typical 
coronal active region is taken to be 2 x 10^ km, indicates that for all 
three transitions the optical depth in the Doppler core is about two. 

This should not affect these results significantly, however, since if we 
assume that at most half of the resonantly-scattered photons are directed 
towards the chromosphere or photosphere (where they would be absorbed), 
then the total emission would be reduced by no more than ~ 407,. For 
these three lines, the values of KG(T) were calculated using van Regemorter's 
values of P(W/kT), Cross ley and Dalgarno's (1965) f-values, and Jordan's 
(1969) tabulations of n(ion) /n(element) . n(Fe)/n(H) was taken to be 
5 x 10‘ 5 , the value used in Tucker and Koren's work. 

There is a good correspondence between F 0 b s and F ca x c > as can be 
seen from the plot of Figure 1. The differences between F Q ^ g and F ca -j_ c 
are mostly within a factor of two or so. This is acceptable in view of 
the uncertainties in the atomic parameters. Discrepancies much larger 
than this (e.g. the lines of 0 VII at 21.60A and 0 VIII at 14.82A, fluxes 
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of which were taken from Evans and Pounds (1968)), seem to be explicable 
in terms of possible errors in the observational data. The 0 VII line, 
for example, occurs in a region where the efficiency of Evans and Pounds 1 
detector was very low, while the 0 VIII line had a very low count rate. 

The three ultraviolet lines due to Fe XV and Fe XVI are rather more intense 
than is calculated from Chambe r s model, not less as would be expected if 
some of the resonantly- scattered photons were directed downwards and ab- 
sorbed by the photosphere. Since the observation of these lines were 
taken on an occasion different from those of the x-ray lines, there is a 
possibility that the enhanced ultraviolet emission arose from a different 
level of solar activity. In conclusion, then, the form of 0(T) given by 
Chambe would seem to describe non- flaring active regions rather closely. 

Some further evidence that the form of 0(T) which Chambe discusses 
is close to the correct form is presented by Walker et al . Their functional 
form has a similar shape to Chambe ! s, and with it they are able to repro- 
duce the observed fluxes of eight lines between 4 and 8 A very well indeed. 
However, this has been done through adjustment of element abundances 
which appears to affect the analysis to a greater extent than the precise 
form of 0(T). Moreover, although the continuum flux at X = 8A is correctly 
predicted by their model, that at the shorter wavelengths is not, sug- 
gesting a greater quantity of emitting material at T > 7 x 10°K than is 
suggested in their model. 

5 . SOME OTHER METHODS FOR THE SOLUTION OF 0(T) 

In this section, we discuss two methods for the derivation of the 
functional form of 0(T) which hitherto have not been tried but which ought 
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to have an advantage over Chambe's method in that a determination of 0(T) 
is made ab initio . Both procedures enable a large number of spectral 
lines to be used with a least-squares analysis, so that there is some 
hope that the effect of poorly known atomic parameters may be reduced to 
some extent. 

(a) 0(T) Taken to be a Power Series in T . 

At first, we shall take 0(T) to be given by 
+N2 

0(T) = S a n T n (Ni, No ^ 0). (9) 

n=-N 1 * W 


With a set of m line fluxes (m greater than or equal to the number of 
terms in the summation), F^(i ~ 1, ..., m), we can solve for 0(T) from 
(9) and (8) using 


n 2 

Fi = K i E a J* G ± CT> T n dT, i - 1, 
n=-N x n £ 


* * 0 * j m. 


A least-squares procedure could be used to solve for the unknown coefficients 
a n ra e xceeds the number of terms in the sum. 

There are two difficulties with using the form (9). The first is 
that terms of the summation with n S 3 tend to produce spurious oscil- 
lations which may not be distinguishable from possible real structure in 
0(T). The second is that any representation of 0(T) involving a sum of 
terms may give negative values for 0(T) in certain temperature ranges. 

Using x-ray data for a flare plasma (the results to appear in Paper II of 
this series), it has been found that both effects should be avoidable by 
having a set of lines in which the temperature ranges of the G(T) functions 
(defined, say, by the range within which G(T) is at least 1/ex its maximum 
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value) overlap over the entire temperature interval (T^> T 2 ) of interest. 
The .results of Walker et al . suggest that, for a non-flaring active region, 
temperatures between T 1 = 1 x 10 6 K and T 2 £ 9 x 10 6 K should be considered, 
while 0S0-5 results to be reported in Paper II indicate that, for a flare, 
Tl = 1 x 10^K, T 2 ~ 50 x 10^K. Spurious oscillations may still occur just 
above the temperature interval (Tj, T 2 ) > these may be considerably 

damped by including in the data sample lines which have zero or near-zero 
fluxes and whose G(T) functions extend well beyond the temperature T 2 • 

A series of solutions for 0(T) in the form of (9) can be obtained 
for various values of N b and N 2 • The goodness of fit in each case may be 
tested by comparing the observed line fluxes, F Qbs , with the value of 
K G(T) 0 CO ^1, 0 / (T) being one of the solutions. The integration 
over (T-p T 2 ) should not include any region where 0 / (T) < 0; any solutions 
having this property should be rejected. If the solutions giving best 
agreement between observed and calculated line fluxes show any marked 
similarity with a well-known functional form (e.g. Chambe's expression 
resembles the form 0(T) = A T^ exp (-j3T)), such an expression could be 
tested by varying its parameters to see if there is any further improve- 
ment in reproducing the observed line fluxes. 

The analytic expression for 0(T) which is finally derived may be 
used to find the flux of continuum emission which may then be compared 
with that observed. The x-ray continuum is made up of principally free- 
free and free-bound emission for temperatures typical of non-flaring 
active regions. At wavelengths of about 20— 60A, two-photon emission from 
H-like and He-like ions may also be an important contributor (e.g. Tucker 
and Koren, 1971). Writing the continuum radiation per unit n e ^ dV/dT at 
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energy E (keV) and temperature T as T|(E,T) photons cm - 2 S ~1 keV~^/lpm'^ K - ^), 
the flux of continuous emission predicted by the 0(T) derived as above is 

oo 

^(E) = J 0CT) T|(E,T) dT photons cm“^ s - ^ keV~k 
o 

This may then be compared with observed values of ^(E) at spectral locations 
relatively free of line emission. 

(b) Gauss-Laguerre Quadrature Method 

In this method, we replace the integral of (8) with a quadrature 
sum. We substitute x ■= cT (c = constant) for the integrating variable, T, 
so that (8) becomes 

K “ 

F = ~ J G(x) 0(x) dx. (10) 

o 

Writing 

f(x) = G(x) 0(x) e x , ( 11 ) 

and replacing the infinite integral of (10) by a Gauss - La guerre quadrature 
sum, we have 

K a 

F * c 2 A k f(x k ) (12) 

k=l 

where the x k s are zeros of the nth Laguerre polynomial and the A k s weight 
factors. Replacing the integral of (10) with the quadrature sum of (12) 
is exact if f(x) is a polynomial in x of degree < 2n - 1. Values for x^ 
and A k are tabulated for n = 1,2, 3,..., 16 by Krylov (1962). Writing (12) 
in terms of temperature, we have for the flux of the ith line of a sample 
of m lines, 


*i “ h. ^ A k Gi (T k ) 0(T k ) e cT k. 
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If m > n, we may then solve for the n unknowns, 0(T^) (k-1 ,2 , . . . ,n) , 
using a least-squares procedure having chosen some suitable value for 
the constant, c. 

The method thus finds discrete values for 0(T). It may not be suit 
able if the function f(x) in ( 11 ) cannot be adequately described by a 
polynomial of degree < 2n - 1. Values of 0(T^) may only be established 
if all the are within the temperature range of the various G(T) 
functions for the m lines* 


6. SUMMARY 

The functional form of 0(T) = n e 2 dV/dT can in principle be derived 
for a coronal active region using intensities of x-ray or ultra- 

violet lines from optically thin plasma. The procedure assumes that the 
electron density, ne = r^CT), is a single-valued function of temperature. 

A large number of lines will give a better determination of 0(T) if all 
the atomic parameters and the observed fluxes are reasonably well-known. 

At the present time, errors may still be significant for the rates of 
electron collisional excitation and cascade contributions, particularly 
for ions more complex than the He-like stage. Errors in values for element 
abundances (arising from the possibility of heavy-element enhancement by 
diffusion) and ion concentrations may also have to be considered. 

Previous attempts to determine 0(T) are examined critically. The 
least-squares method of Batstone et al . (1970) is based on very low 
photon count rates, a fact which appears to explain why some values of 
jT+AT 0 ( T ) ^ are ne g a tive when no constraints are imposed on the analysis. 
The results of the re -analysis are summarized in Table I. Chambe's (1971) 
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method was to take some arbitrary functional form for 0(T) and to adjust 
its shape until a good fit with observations is achieved. It is not 
shown, however , that the derived 0(T) actually reproduces the observed 
line fluxes. The atomic parameters used in the analysis, moreover, were 
in many cases only rough approximations. A re-analysis is made in which 
the x-ray line fluxes are calculated from (8) using the 0(T) as given by 
Chambe and more accurate atomic data. The agreement between the observed 
fluxes and those calculated as a result is mostly to within a factor of 
two, so that we must conclude that Chambers form for 0(T) is able to de- 
scribe the x-ray line emission rather well. 

Two ab initio methods are developed in which 0(T) may be derived 
from x-ray line fluxes. In the first, 0(T) is expressed as a power series 
in T, with coefficients to be determined by least-squares methods if the 
number of lines is sufficiently great. That solution (0'(T)) which is 
able best to reproduce the observed line fluxes should be selected. If 
this 0'(T) resembles some familiar analytic form, then this form should 
be tested to see whether any improvement in reproducing observed fluxes 
can be obtained. 

The second method approximates the integral of (8) with a Gauss- 
Laguerre quadrature sum with n terms. This procedure is exact if the 
right side of (11) can be expressed as a polynomial of degree < 2n - 1. 
With a set of m (> n) x-ray lines, solutions for the n unknown values of 
0(T^.), k - 1, n, may be obtained by a least squares procedure. 
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FIGURE CAPTION 


Plot of observed flux (erg cnT'^ s“*-) v£. flux calculated from Chamber 
0(T) function. Points are annotated with ion and line wavelength in A- 



OBSERVED LINE FLUX - erg 





